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Abstract 

We calculate the intensity of the polariton mediated inelastic light scattering in semiconductor 
microcavities. We treat the exciton-photon coupling nonperturbatively and incorporate lifetime 
effects in both excitons and photons, and a coupling of the photons to the electron-hole continuum. 
Taking the matrix elements as fitting parameters, the results are in excellent agreement with 
measured Raman intensities due to optical phonons resonant with the upper polariton branches in 
II- VI microcavities with embedded CdTe quantum wells. 

PACS numbers: 71.36.+C, 78.30.Fs, 78.30.-j 
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Planar semiconductors microcavities (MC's) have attracted much attention in the last 
decade as they provide a novel means to study, enhance and control the interaction between 
light and matter . When the MC mode (cavity-photon) is tuned in near 

resonance with the embedded quantum-well (QW) exciton transitions, and the damping 
processes involved are weak in comparison to the photon-matter interaction, the eigenstates 
of the system become mixed exciton-photon states, cavity-polaritons, which are in part light 
and in part matter bosonic quasi-particles l|, l2|, yj. Examples of interesting new physics are 
the recent evidence of a Bose- Einstein condensation of polaritons in CdTe MC's js, ^, and 
the construction of devices which increase the interaction of sound and light, opening the 
possibility of realizing a coherent monochromatic source of acoustic phonons Tj]. 

Raman scattering due to longitudinal optical (LO) phonons, being a coherent process is 
intrinsically connected with the cavity-polariton. The physics of strongly coupled photons 
and excitons, the polariton- pho non interaction, and the polariton-external-photon coupling 



are clearly displayed 



0, y, m 
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12l |. In particular resonant Raman scattering (RRS) 



experiments, in which the wave-length of the incoming radiation is tuned in a way such that, 
after the emission of a LO phonon, the energy of the outgoing radiation coincides with that of 
the cavity polariton branches, have proven to be suited to sense the dynamics of the coupled 
modes, and to obtain information about the dephasing of the resonant polaritonic state 
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12l |. Unfortunately, due to the large remaining luminesence, RRS experiments 



in resonance with the lower polariton branch have not yet been achieved in intrinsic II- VI 
MC's 131]. Therefore all reported experiments in these kind of MC's, with embedded CdTe 
QW's, consider the case of the scattered photons in outgoing resonance with the upper 
polariton branch (for two branch-systems: coupling of one exciton mode and the cavity 
photon) or with the middle polariton branch (for tree branch-systems: coupling of two 



^0 



14| . The measured intensities in these two 



exciton modes and the cavity photon) [9|, 
systems were analyzed on the basis of a model in which one (two) exciton states |e) are 
mixed with a photon state |/) with the same in-plane wave vector k, leading to a 2x2 (3x3) 
matrix. Essentially, the Raman intensity is proportional to 



I ^ T,TMPAH'\P. 



(1) 



where Tj describe the probability of conversion of an incident photon into the polariton 
state |Pi), Tg has an analogous meaning for the scattered polariton \Ps) and the outgoing 
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photon \fs), and H' is the interaction between electrons and the LO phonons. At the 
conditions of resonance with the outgoing polariton, T, is very weakly dependent on laser 
energy or detuning (difference between photon and exciton energies), while Tg is proportional 
to the photon strength of the scattered polariton \ {fs\Ps)\'^- Similarly, in the simplest case 
of only one exciton (2x2 matrix) one expects that the matrix element entering Eq. ([T]) is 
proportional to the exciton part of the scattered polariton | (Pj|if'|Ps) P ~ K^sl-^s)!^- Thus 
if the wave function is \Ps) = a\fs) + Pl^s), one has 

I-\{fs\Ps)\'\{es\Ps)\' = \am'. (2) 

This model predicts a Raman intensity which is symmetric with detuning and is maximum 
at zero detuning. In other words, the intensity is maximum for detunings such that the 
scattered polariton is more easily coupled to the external photons, but at the same time 
when the polariton is more easily coupled to the optical phonons, which requires a large 
matter (exciton) component of the polariton. This result is in qualitative agreement with 



experiment 10]. However, for positive detuning the experimental results fall bellow the 
values predicted by Eq. ([2]) (see Fig. [1]). This is ascribed to the effects of the electron- hole 
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continuum above the exciton energy, which are not included in the model 

For another sample in which two excitons are involved, the above analysis can be ex- 
tended straightforwardly and the intensity depends on the amplitudes of a 3x3 matrix and 



matrix element of H' involving both excitons lOj. However, comparison with experiments at 



resonance with the middle polariton branch, shows a poorer agreement than in the previous 
case (Fig. 6 of Ref. Q). In addition, a loss of coherence between the scattering of both 
excitons with the LO phonon was assumed, which is hard to justify. Some improvement has 
been obtained recently when damping effects are introduced phenomenologically as imagi- 
nary parts of the photon and exciton energies, but still a complete loss of coherence resulted 
from the fit Q,[l5|. 

In this paper we include the states of the electron-hole continuum, and the damping 
effects in a more rigorous way. Using some matrix elements as free parameters, we can 
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describe accurately the Raman intensities for both samples studied in Ref. 

In the experiments with CdTe QW's inside II- VI MC's, the light incides perpendicular to 
the {x, y) plane of the QW's and is collected in the same direction z. Therefore the in-plane 
wave vector K = 0, the polarization of the electric field should lie in the {x, y) plane, and 
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the excitons which couple with the hght should have the same symmetry as the electric 



field (one of the two Fs states of heavy hole excitons 16|). Thus, to lighten the notation 
we suppress wave vector and polarization indices. The basic ingredients of the theory are 
two or three strongly coupled boson modes, one for the light MC eigenmode with boson 
creation operator /"'", another one for the Is exciton (ei) and if is necessary, the 2s exciton 
(62) is also included. We assume that each of these boson states mixes with a continuum of 
bosonic excitations which broadens its spectral density. In addition, we include the electron- 
hole continuum above the exciton states, described by bosonic operators cj^, Ck, where k is 
the difference between electron and hole momentum in the (x, y) plane (the sum is K = 
because it is conserved). 
The Hamiltonian reads: 

H = Efff + E^eh^ + Y.^v4f + H.c.) + ^ e.rjr, + 5^( Vj/ + H-c.) 

i i p p 

+ e^q4^^q + 5^(^iq4e, + H.c.) + ^ CkcU + 5^(V^.4/ + H.c). (3) 

iq iq k k 

The first three terms describe the strong coupling between the MC photon and the exciton(s) 
already included in previous approaches 10|, ll5|. The fourth and fifth terms describe a 



continuum of radiative modes and its coupling to the MC light eigenmode. Their main effect 
is to broaden the spectral density of the latter even in the absence of light-matter interaction. 
The following two terms have a similar effect for the exciton mode(s). The detailed structure 
of the states described by the dl^ operators is not important in what follows. They might 
describe combined excitations due to scattering with acoustical phonons. The last two terms 
correspond to the energy of the electron hole excitations and their coupling to the MC light 
mode. 

In Eq. (IH]) we are making the usual approximation of neglecting the internal fermionic 
structure of the excitons and electron-hole operators and taking them as free bosons. This 
is an excellent approximation for the conditions of the experiment. We also neglect terms 
which do not conserve the number of bosons. Their effect is small for the energies of 
interest jl?!]. These approximations allow us formally to diagonalize the Hamiltonian by a 
Bogoliubov transformation. The diagonalized Hamiltonian has the form H = 'YliyE^plpy, 
where the boson operators pi correspond to generalized polariton operators and are linear 
combinations of all creation operators entering Eq. (I3l). Denoting the latter for brevity as fej. 
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then pI = J2j ^ujb\- In practice, instead of calculating the Ey and Ayj, it is more convenient 
to work with retarded Green's functions Gji{uj) = {{bj] h\))^ and their equations of motion 

u{{h,-,h\))^ = 5,i + {{[h,,Hlh\))^. (4) 

As we show below, the RRS intensity can be expressed in terms of spectral densities derived 
from these Green's functions, which in turn can be calculated using Eq. (jl]). 

Using Fermi's golden rule, the probability per unit time for a transition from a polariton 
state \i) = pIi\0) to states \s) = pla^O), where creates a LO phonon is 

W ='^\{i\H'\s)\'^p{u;), where p(cu) = y5(cu-E^) = --ImG^^(w + zO+) (5) 

is the density of final states. As argued above, we neglect the dependence of Tj on frequency 
and take Tg = |v4j,jp, the weight of photons in the scattered eigenstates. In addition, H' 
should be proportional to the matter (exciton) part of the scattered polariton states. Then, 
the Raman intensity is proportional to: 

WTs^I= \A,ei + aA,e2\^\A,f\^p{uj). (6) 

Here we are neglecting the contribution of the electron-hole continuum to H', and a is the 
ratio of matrix elements of the exciton-LO phonon interaction between 2s and Is excitons. 
If the 2s excitons are unimportant, a = 0. 
Using Eqs. (jlj) it can be shown that 

pAuj) = -^[G,i{uj + iO+) - Gji{uj - iO-)] = A,,A,ip{uj). (7) 

ZTT 

From here and J2j l^i^jP = 1, it follows that p{uj) = J2jPjji^)- Replacing in Eq. we 
obtain: 

7/ N P//[Pel,el + |aPPe2,e2 + 2Re(ape2,el)] /qx 

'e;^) ■ 

In practice, when u is chosen such that the resonance condition for the outgoing polariton 
is fulfilled, we can neglect the contribution of the continuum states in the denominator of 
Eq. ([8]). In particular, if the contribution of the 2s exciton can be neglected (as in sample 



A of Ref. [lO) 



liuj) = PfAfP^^A-^) (9) 
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The Green's functions are calculated from the equations of motion (jlj). In the final expres- 
sions, the continuum states enter through the following sums: 

sjuj) = y — — , s,icu) = y — — , s'{uj) = y — — (lo) 

p ^ q ^ k 

For the first two we assume that the results are imaginary constants that we take as param- 
eters: 

Sf{uj) = —iSf, Sj{uj) = —i6j (11) 

This is the result expected for constant density of states and matrix elements. Our results 
seem to indicate that this assumption is valid for the upper and middle polariton branches. 
For the lower branch at small k it has been shown that the line width due to the interaction 
of polaritons with acoustic phonons depends on detuning Ef — Ei and k\\, being smaller for 
small wave vector Q ■ 

The electron-hole continuum begins at the energy of the gap and corresponds to vertical 
transitions in which the light promotes a valence electron with 2D wave vector k to the 
conduction band with the same wave vector. In the effective-mass approximation, the energy 
is quadratic with k and this leads to a constant density of states beginning at the gap. The 
matrix element is proportional to Mk = (kt,|p£;|kc), where pe is the momentum operator 
in the direction of the electric field, and |k.t,), |kc) are the wave functions for valence and 
conduction electrons with wave vector k. Taking these wave functions as plane waves, one 
has Mk ^ kE, the wave vector in the direction of the electric field. Then |VkP ~ k"^- Adding 
the contributions of all directions of k one has |VkP ~ |kp ~ (linear with energy for small 
energy). This leads to 

S'f = R{uj) -iA{uj - B)e{uj - B), (12) 

where B is the bottom of the electron-hole continuum (the energy of the semiconductor gap) 
and A is a dimensionless parameter that controls the magnitude of the interaction. The real 
part R{uj) can be absorbed in a renormalization of the photon energy and is unimportant 
in what follows. The imaginary part is a correction to the photon width for energies above 
the bottom of the continuum. 



Using the theory outlined above, we ca^ 



the samples A and B measured in Refs. 
results. 
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culated the intensity of RRS corresponding to 



14 and compared them with the experimental 
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Sample A corresponds to the simplest case. Two polariton branches are seen and therefore 
only the Is exciton plays a significant role. The binding energy of this exciton B — Ei is not 
well known. The Rabi splitting 2Vi = 19 meV. The width of the Raman scan as a function 
of frequency for zero detuning is of the order of w = 0.1 meV (see Fig. 3 of Ref. llOl ). This 
implies the relation = 6'j + 6l in our theory. In any case the results are weakly sensitive 
to w. Therefore, we have three free parameters in our theory in addition to a multiplicative 
constant: B — Ei, the ratio of widths Sf/6i and the slope A. 




FIG. 1: Raman intensity as a function of detuning for sample A. Solid squares: experimental 
results Q- Solid line: theory (Eq. ^ foi B - Ei = 14 meV, 6f = 6i= 0.1 meV, and A = 0.031. 
Dashed line: result for a 2x2 matrix (Eq. [2]). 

The comparison between the experimental and the theoretical intensities is shown in 
Fig. [T]for a set of parameters that lead to a close agreement with experiment. The condition 
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of resonance is established choosing the energy uj for which the intensity given by Eq. ([9]) has 
its second relative maximum (corresponding to the upper polariton branch). The dashed line 
corresponds to the case in which only the first three terms (with i = 1) in the Hamiltonian, 
Eq. ([3]) are included. In this case, the intensity is given in terms of the solution of a 2x2 



matrix [Eq. ([2])] and was used in Ref. [id to interpret the data. This simple expression gives 
a Raman intensity which is an even function of detuning. When the full model is considered, 
the Raman intensity falls more rapidly for large detuning Ef — Ei as a consequence of the 
hybridization of the photon with the electron-hole continuum. When the energy of the 
polariton increases beyond B entering the electron-hole continuum (corresponding to the 
kink in Fig. 1), the Raman peak broadens and loses intensity. The kink can be smoothed 
if the effect of the infinite excitonic levels below the continuum is included in the model 
(leading to an 5*^ with continuous first derivative), but this is beyond the scope of this 
work. If the ratio 5f/5i is enlarged, the Raman intensity increases for negative detuning 
with respect to its value for positive detuning. 



In the experiments with sample B three polariton branches are observed iJ] and the 
2s exciton plays a role. Experimentally, it is known that the binding energy for the two 
excitons are B — Ei =17 meV and B — E2 =2 meV. From the observed Rabi splitting one 
has 2Vi = 13 meV, 2V2 = 2.5 meV. In comparison with the previous case, we have the 
additional parameter a (the ratio of exciton-LO phonon matrix elements). In addition, to 
be able to describe well the intensity for low energies of the middle polariton (left part of 
the curve shown in Fig. [2]), we need to assume a small linear dependence of Ei with the 
position of the incident laser spot in the sample. This dependence is also inferred form the 
observed luminescence spectrum [l^. In our model, this corresponds to a dependence of Ei 
with Ef-. 

El = E° + z{Ef - El) (13) 

In Fig. [2] we show the intensity at the second maximum of I {00) (corresponding to the 
middle polariton branch) as a function of the energy of this maximum. We also show in the 



figure experimental results taken at 
(4.5 K) than those reported in Ref. 



ower laser excitation and a slightly higher temperature 



10. 



The slope which better describes the data is z = 0.14. This value is close io z = 0.155 
which was obtained from a fit of the maxima of luminescence spectrum of the lower and 
middle polaritons. We have taken the same value for A as in Fig. [TJ The agreement between 
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theory and experiment is remarkable. As for sample A, 
fit are reasonable in comparison with calculated values 



;he values of Si that result from the 
3. 




FIG. 2: Raman intensity as a function of the middle polariton energy. Solid squares: experimental 



results 



Solid line: theory for 6f = 0.2 meV, 5i = 0.1 meV, 62 = 0.12 meV, a = —0.45 and 



z = 0.14. A is the same as in Fig. [T] . 

In summary, we have proposed a theory to calculate Raman intensity for excitation 
of longitudinal optical phonons in microcavities, in which different matrix elements are 
incorporated as parameters of the model. The most important advance in comparison with 
previous simplified theories [lo|, [l^ is the inclusion of the strong coupling of the electron- 
hole continuum with the microcavity photon. Inclusion of this coupling is essential when the 
energy of the polariton is near the bottom of the conduction band (at the right of Fig. 1). 
We also have included the effects of damping of excitons and photons, coupling them with 



9 



a continuum of bosonic excitations. Simpler approaches have included the spectral widths 
5f and 5i of photons and excitons as imaginary parts of the respective energies, leading to 
non-hermitian matrices. 

Taking some of the parameters of the model as free {6f/6i, B — Ei and A for sample A, 
5f, 6i, z and a for sample B), we obtain excellent fits of the observed Raman intensities. The 
resulting values of the parameters agree with previous estimates, if they are available. We 
are not aware of previous estimates for A and a. As an important improvement to previous 
approaches Q, Q for the case of sample B, we do not have to assume a partial loss of 
coherence between Is and 2s excitons in their scattering with the LO phonon. 

Further progress in the understanding of the interaction of excitons with light and 
phonons requires microscopic calculations of the parameters 6f,6i,A and a, and the effects 
of the temperature on them. However, taking into account the difficulties in calculating 
these parameters accurately, the present results are encouraging and suggest that the main 
physical ingredients are included in our model. 
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CONICET and PICT 03-13829 of ANPCyT. 
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